program run 
use cdata
use technical
use df_dt
use timestep
implicit none
logical :: lexist
integer :: ierr
real,allocatable,dimension(:) :: y
!<define namelists > 
namelist /common_pars/ &
    ndim,model
!<First read in the input files >
  inquire(FILE="ode.in", EXIST=lexist)
if (.not. lexist) call fatal_error('ode.in','input file does not exist') 
!--------------------------------------------------!
open(unit=1,FILE='ode.in',FORM='formatted',STATUS='old')
read(1,NML=common_pars, IOSTAT=ierr)
!rewind(1)
if(ierr/=0) call fatal_error('common_pars','problem with reading common_pars')
!<Allocate the arrays needed>
allocate(y(ndim))
!<initialize >
call ini_rhs(model,1)
!<close the file>
close(1)
!---------------------------------------
deallocate(y)
!---------------------------------------
end program
